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Abstract 

We advance scale-invariance arguments for systems that are governed (or approximated) by a q— Gaussian 
distribution, i.e., a power law distribution with exponent Q = 1/(1 — q); q 6 TSL. The ensuing line of 
reasoning is then compared with that applying for Gaussian distributions, with emphasis on dimensional 
considerations. In particular, a Gaussian system may be part of a larger system that is not Gaussian, but, 
if the larger system is spherically invariant, then it is necessarily Gaussian again. We show that this result 
extends to q-Gaussian systems via elliptic invariance. The problem of estimating the appropriate value for 
q is revisited. A kinetic application is also provided. 
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I. INTRODUCTION AND BACKGROUND MATERIAL 



Homogeneous power-laws, such as Newton's universal law of gravitational attraction, for instance, 
abound in Nature. They are, by definition, self-similar and thus true in all scales. Systems sta- 
tistically described by power-law probability distributions are rather ubiquitous [1] and thus of 



perennial interest yfl. In this report we wish to give careful scale-invariance consideration to 
systems that are governed (or described) by a special kind of power-law probability distribution 
functions (PDF), namely, the q— Gaussian function. Consider a system S described by a vector 
X with n components. We say that X is q— Gaussian distributed if its probability distribution 
function writes as described by (|L?1)-(|L6l) below. 

It is well-known that for such S— systems one can appeal to Jaynes' maximum entropy principle 
(MaxEnt) yO under a covariance constraint with a generalized (or q— ) information measure 

H q 0) = (1 - q)- 1 I dx [f{x) - q G K, (1.1) 

n 

as the protagonist [|4J]. This measure has been found to be useful in extracting information per- 
taining to systems that are characterized by either (1) fractal nature, (2) long-memory, or (3) long- 
range interactions [5|]. Employing H q for other types of system has generated controversy [60 
which is of no relevance to our present purposes. We will focus attention upon properties of Gaus- 
sian systems that remain valid for their q-Gaussian counterparts as well (as q becomes different 
from unity), with emphasis on the dimensional properties of both kinds of systems. It is well- 
known that if a system "is Gaussian", any part (sub-system) of it is still Gaussian. This property 
holds for q— Gaussian systems as well, as proved in [7]. A more interesting result is the inverse 
phenomenon: a Gaussian system may be part of a larger system that is not Gaussian. However, 
if the larger system is spherically invariant, then it is necessarily Gaussian again. Surprisingly 
enough, this "inverse property" has gone largely ignored in the statistical literature. In this work 
we will not only provide a simple proof for it but we will show that it can be extended to q-Gaussian 
distributions as well. These results can be given a physical interpretation within the framework of 
the estimation of the parameter q of a given system [50. We will prove that, if spherical symmetry 
prevails, such estimation can be performed using only a restricted, observable part of the system 
and that the overall parameter q for the entire system can be retrieved provided the dimension of 
the system is known. We begin our considerations by introducing the two basic notions, namely 
spherical symmetry and g— Gaussian systems. We will also apply our results to a simple case of 



kinetic theory via the Beck-Cohen superstatistics theory 



ffiHHHHQ. 



A. Spherical symmetry 



A really momentous symmetry is that of invariance against rotations. It is found in the fun- 
damental laws of nature and constitutes one of the most powerful principles in elucidating the 
structure of individual atoms, complicated molecules, and entire crystals. Also, it characterizes 
the shape of many systems. We can cite self-gravitating systems like stars and planets, that have 
quasi- spherical shape if their mass is large enough. Also, many atomic nuclei are spherical, and 
many molecules as well, etc. Conservation of angular momentum, a very frequent occurrence, is a 



result of the isotropy of space itself 111 511 . We discuss now some properties of spherical probability 
distributions. 

The characteristic function associated with a random vector X E W 1 is 

<p x (U) = Ee lUtx ; U el" (1.2) 

Under the hypothesis, discussed for instance in the textbook [16, XV.3], that cp x G £i(M n ), there 
is a one-to-one relation between ip x and the probability density function f x of X. The random 
vector X is said to have a spherical distribution if its characteristic function ip x satisfies 

<Px{U) = <KW\\) (1.3) 

for some scalar function <p : IR + — > K which is then called the characteristic generator of the 
spherical distribution. We will write X ~ S n ((f>) in this case. It is well-known that an equivalent 
definition for a spherical random vector X is 

X ~ AX; V A orthogonal 

where ~ denotes equality in distribution. 

Spherical random vectors have, as it is well-known, the following properties: 

1 . All marginal distributions of a spherical distributed random vector are spherical. 

2. All marginal characteristic functions have the same characteristic generator. 

3. If X ~ S n (<f>) then 

X ~ rT n (1.4) 

where T n is a random vector distributed uniformly on the unit sphere surface in W 1 and r is 
a positive random variable independent of T n . 



Let us remark that a spherically distributed random vector does not necessarily possess a density. 

A generalization of the concept of spherical distribution is given by elliptical distributions, to 
which the multi-normal distribution belongs. Elliptical distributions have recently gained a lot of 
attention in financial mathematics, being of use particularly in risk management. A n— random 
vector Y is said to have an elliptical distribution with so-called characteristic matrix Cy{n x n) 
if Y ~ AX, where X ~ <S n (0) and A is a n x n deterministic matrix such that A 1 A = Cy and 
rank(Cy) = n. We shall write Y ~ £ n (Cy, 0). 



B. q— Gaussian systems 



An n— components vector X is q— Gaussian distributed if its PDF writes as follows |2|] 



in the case 1 < q < 

^ n+2 



f x (X) = A q (l + X t A- 1 X) - 



matrix A being related to the covariance matrix K = EXX 1 in the fashion fl 1 Vfl 



(1.5) 



A = (m - 2) K. 



where the number of degrees of freedom is defined as QJ/Ztl 

2 



m 



q-1 



— n. 



Moreover, the partition function Z q = 1 /A q reads QJ/ZD 



r * 



and the characteristic function is 



(1.6) 



(1.7) 



<Px(U) 



2 x "t 



Z 2 Kni(z) 



with z = VWAU and K is the modified Bessel function of the second kind. 



(1.8) 



in the case q < 1 



f x (X) = A q (l-X t ^- l X)X 



(1.9) 
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with matrix £ = dK and parameter d defined as d = + n. In this case, the partition 
function is 

r (fEf) kE|v a 



r f -i- 2 

H T 2 



and the characteristic function is 



d 1 d, Ji-l\ z ) 

<Px (U) = 2^ 1 r(-)s- rrr - (i.io) 

Z2~ 



where z = yTPTJJ and J is the Bessel function of the first kind. 
We begin in the next Section to advance our present results. 

H. SIZE BEHAVIOR OF THE q —PARAMETER 

Our first result revolves around the behavior of the non-extensivity parameter q as a function 
of the dimension of the system and is embodied in the following theorem, the proof of which is 
given in the Appendix. 

Theorem 1 Assume that a system X n G W 1 follows a q— Gaussian distribution with parameter q n ; 
then with 1 < k < n — 1, any k— dimensional subsystem X k = [xi, . . . , Xk\ of X n is q— Gaussian 
distributed with parameter 

* =1 -2 + ( B -t)(l- fc ) <IU) 
Reciprocally, assume that a n— dimensional and spherical system X n contains a 
k— dimensional subsystem X k that follows a q— Gaussian distribution with parameter q k ; then 
system X n is itself q— Gaussian distributed with parameter q n defined as in (III.U 

Recall that for any n— dimensional orthogonal transformation A, there exists an orthogonal 
decomposition of M n as 

R" = E © F © Gt © • • • © G k 
into stable subspaces, the restriction of A to each subspace being 

• the identity transformation for subspace E 

• minus the identity transformation for subspace F 



• a two-dimensional planar rotation of angle 6 k for subspace G k 

Moreover, n— dimensional q— Gaussian distributions with parameter q n arise in statistical physics 
as the canonical distributions of systems with maximal q— entropy H g of order q n and fixed covari- 
ance matrix. The result of theorem CO can be paraphrased in this way: if a system of dimension k 
has maximum g— entropy of order q k , and is also part of a spherical system of dimension n > k, 
then the whole system maximizes the q— entropy with parameter q n related to qk as in (|II.1I) . No- 
tice that in (III. II ), as q k — > 1 then q n — > 1 and we deduce that if a spherical system has a Gaussian 
part, then it is necessarily Gaussian as well. 

We note also that q n < 1 implies q k < 1 while q n > 1 implies q k > 1. This result is natural 
since cases q > 1 and g < 1 correspond to two different t ype s of distributions: according to Beck 
and Cohen's superstatistics principle [|8|,|9|, 12LH3LHJ: 



• g > 1 that corresponds to a Gaussian system subjected to fluctuations that are independent 
of the state of the system. 

• q < 1 that corresponds again to a fluctuating Gaussian system for which the amplitude of 
the fluctuations depends on the system's state. 

Alternatively, these two cases can be characterized as follows. Our system is here described by 
the random vector X: if q > 1, then X has unbounded support, contrarily to the bounded support 
associated to the case q < 1 II 1711 . 

Thus, we obtain the following rather natural result: reduction or enlargement of a q— Gaussian 
system does not change its superstatistical nature. 

IE. A SECOND RESULT: MULTI-COMPONENT SYSTEMS 
A. Average behaviour 

It may happen that measuring the behavior of only one component X\ of a large system lei™ 
is not physically feasible, and one has to content oneself with measuring instead the behavior of a 



superposition of contributions from (or average of) several components lll8ll . in the form 



< x >=^2a iXi (DTI) 



1=1 
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where the deterministic coefficients a; 6 1 characterize the measurement device. The following 
theorem (see Ref. JljJ]) allows to give a special characterization of this average value in the case 
of spherical systems. 

Theorem 2 [19, Th. 2.4] If X e W n is spherically distributed and A = (a 1; . . . , a n )* is a deter- 
ministic vector then 

n 

< X >= CLiXi 
i=l 

is distributed as \\A\\xi, where \\A\\ is the Euclidean norm of A. 

In the next subsection we extract rather interesting physical conclusions from this theorem. 



B. Application to the estimation of q 

An important problem in non-extensive statistics is the estimation of the non-extensivity pa- 
rameter q n associated to an n— dimensional system that follows a q— Gaussian distribution [5]. We 
provide here some hints about a possible estimation strategy in the case q > 1, assuming that we 
have access to averaged measures of the system of the type QUI .lb . 

Assuming q n > 1, then if X n follows distribution (OT) . the averaged measure < X > is 
distributed as 

A a , . X * m+l 



f<x>(x) = tttUI + 



2 



with m = 2/(qi — 1) — 1 and A = A^i. As a consequence, a possible estimation strategy of 
parameter q n follows the three following steps: 

1 . since the "tail-behavior" of the distribution of /<x> is 

f <x> (x) ~ x~^ +1 \ 

(where ~ means here asymptotic equivalence), parameter m can be estimated as the Levy 
exponent of the distribution of the average measure of the system I2O I 



2. the non-extensivity parameter q 1 of < X > can be computed using (|L7I) as 

m + 3 
m + l 
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3. the non-extensivity parameter q n of the n— dimensional system X n can in turn be deduced 
using (III. II) as 

_ 2-(n + l)(l- gl ) 
qn 2-(n-l)(l- ft ) 

As a new result we find that if the dimension n of the system is known, its non-extensivity 
parameter q n > 1 can be evaluated from any measurement of the type ( ini.lt . 



C. A kinetic application 

1. Theoretical framework 

Another application of the latter result can be provided in the context of the kinematics of 
collision events. We envision a scenario in which attention is focused on the particles of a system 
interacting with a heat bath (a fundamental problem in thermodynamics). An elastic collision 
between (i) a system's particle with momentum P, mass M, velocity V, and energy E and (ii) a 
particle from the heat bath with momentum p, mass m, velocity v, and energy e, verifies J21I 



E + e = E + e 
P + p = P + p 

where "hats" refer to quantities after the collision. In the non-relativistic case, these quantities 
write 

p = mv 

where momenta are 3— dimensional quantities. These equations can be solved as 

Assuming that P and p are independent random variables, we look for stationary distributions 
for p and P, that is, for probability density functions f p and f P such that if p ~ f p and P ~ f P 



P 
E 



MV, 
\\p\\ 2 



2M ' 
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then after the collision, p ~ f p and P ~ f P . An obvious pair of stationary solutions is given 112 lh 

by the independent Maxwell solutions 

1 / IIPII 2 \ 1 / \\pP 

fp( p ) = ~ — — -— ^exp -— — — , Up) = - — — — exp 



(27rMk B Tf 2 * \ 2Mk B Tj ' (27rmk B Tf 2 * \ 2mk B T 



(III.2) 



2. The correlated scenario 



Suppose however that the assumption of independence between momenta p and P does not 
hold. Such is the case, for example, when the corresponding particles are subject to the same 
fluctuations (an interpretation for this scenario is provided in the following Subsection). In this 
instance we look for a stationary joint distribution for p and P, i.e., for a probability density 
function f Pt p such that if (p, P) ~ / Pi p then, after the collision, (p, P) ~ f Pt p again. We note that 
this in turn implies p ~ f p and P ~ f P . 

We need here an extension of Theorem [2] as given below, the proof of which can be found, for 
example, in J22 I. 

Proposition 3 If X ~ S n (Cx, <P) and A is a full-rank (n x n) matrix then Y = AX ~ £ n (Cy, 0) 
with 

C Y = AC x A l . 



Now assume that 

C x 



ml 3 3 
3 Mh 



A 



(m - M)I 3 (2m)/ 3 
(2M)I 3 (M - m)h 



m + M 

where O3 denotes the (3 x 3) null matrix and J 3 the (3 x 3) identity matrix. Then 

C Y = AC x A l = C x - 
We are now in a position to deduce the following result: 



Theorem 4 If I I ~ £q(Cx, <P) with characteristic matrix Cx as above, then the momenta 
P 



vector after the collision 



V 
P 



£&{Cx, 4>). As a consequence, any elliptical joint distribution 



with characteristic matrix Cx is stationary. In particular, p and p have the same distribution, as 
well as P and P. 
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3. Superstatistics at work 



A more physical interpretation can be given to the preceding result, using the notion of super- 
statistics |8J. We know from 112 ill that a pair of independent Gaussian momenta are stationary for 
the collision process. Now, 

• if p ~ A/^mA^T) (the 3— dimensional Gaussian distribution with covariance matrix 
mksTIz) - and 

• P ~ J\f 3 (Mk B T), then 

• p ~ Af 3 (mk B T) and P ~ N 3 (Mk B T). 
Since 

r{ p,f) y M + m )p+ \ M + m) ' 

choosing any (dimensionless) random variable a independent of both p and P, and defining the 
new quantities q = ap, Q = aP, and Q = aP, we deduce that 

Q = ( T7T — I q + \ WjL — ) ^' 

\M+mJ \M + mJ 

so that, obviously, the pair (g, Q) is another couple of momenta whose distribution is stationary. 
Obviously, variables q = ap and Q = aP are not independent ones, since they share the same 
random factor a (unless a is almost surely a constant, which reduces to the Gaussian case). 

As a special case, if a follows an inverse chi-distribution with m degrees of freedom, one 

immediately finds [8J] that the random vector X = | ^ | follows a Tsallis-distribution 

fx (X) =A q (l+ X t A~ 1 X) ^ , (m.3) 



with non-extensivity index q > 1 related to parameter m as in formula (11771) with n — 6. We remark 
that the random variable a can be interpreted, in such a context, as representing temperature's 
fluctuations, as shown by Beck and Cohen Q80 - Thus, the presence of temperature fluctuations 
indicates that the momenta of the incoming colliding particles are correlated. Conversely, if they 
are correlated, then temperature fluctuations ensue. This scenario is a feasible one if the heat bath 
is a finite one, which, in turn, establishes a natural connection with an old result of Plastino and 
Plastino |23|]. 
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IV. CONCLUSIONS 



In this work we have considered physical applications of a largely ignored result of the statisti- 
cal literature: a Gaussian system may be part of a larger system that is not Gaussian. However, if 
this larger system is spherically invariant, then it is necessarily Gaussian again. 

We have provided a simple proof for it and we have shown that it can be extended to 
q— Gaussian distributions as well. Our results have been given a physical interpretation within 
the framework of the problem of estimation of the q-Gaussian parameter q. Also, we applied them 
to a simple instance of kinetic theory involving Beck and Cohen superstatistics |8|]. 

V. APPENDIX: PROOFS 

A. Proof of theorem Q] 

We give here a simple proof of theorem CQ), the principle of which has been kindly suggested 
to us by Pr. Wlodek Bryc and Pr. Jacek Wesolowski. Assuming first that X k = [x%, . . . , x k Y is 
q— Gaussian with parameter q k > 1, we deduce that, with U k = [u\, . . . , UkY, 

VXk (U k ) = <K\\U h \\) = ^\\U k \\fK f (\\U k \\) 

so that 

0(«) = ^NI**?(H) 

and 

~Y Hi 

<PX n {Un) = <P(\\Un\\) = p^ll^ll^dll/nll) 

and X n is q— Gaussian with dimension n and m degrees of freedom, and thus has non-extensivity 
parameter q n such that 

2 t 2 
m = k — n. 

qk - 1 q n - 1 

The same result applies with q k < 1 by considering characteristic function defined as in (11.101) . 

B. An alternate proof of theorem [J 



We provide here an alternate proof based on stochastic representations, as first used in [24], 



extending their result to the case q > 1. Assume that X k e M. k is q— Gaussian distributed with 

11 



parameter q^. a stochastic representation of Xj. is (see Ill9h ) 



-*fc — 

Xd 



where Xk and Xd are independent and chi distributed random variables, where d 
Z k is uniform on the sphere. We know moreover that if X n = rZ n then 



9fc-l 



k and 



where ~ (3 k n^k . We deduce that 



or 



But by Luckacs theorem [25]: 



7 ry Xk 

rdiZ k = — 
Xd 



r 2 d\ = \ 

xl 



we deduce that 



Xn Xk Xk 

xl xl + xl-k xl 



Xn 

Xd 



and that X n = rZ n is q— Gaussian distributed with parameter q n = 1 



2(gfc-l) 



d+n 1 2+(n-k)(q k -l) ' 
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